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Abstract 

Base station (BS) cooperative transmission can improve the spectrum efficiency of cellular systems, 
whereas using which the channels will become asymmetry. In this paper, we study the impact of 
the asymmetry on the performance of channel estimation and precoding in downlink BS cooperative 
multiple-antenna multiple-carrier systems. We first present three linear estimators which jointly estimate 
the channel coefficients from users in different cells with minimum mean square error, robust design 
and least square criterion, and then study the impact of uplink channel asymmetry on their performance. 
It is shown that when the large scale channel information is exploited for channel estimation, using 
non-orthogonal training sequences among users in different cells leads to minor performance loss. Next, 
we analyze the impact of downlink channel asymmetry on the performance of precoding with channel 
estimation errors. Our analysis shows that although the estimation errors of weak cross links are large, 
the resulting rate loss is minor because their contributions are weighted by the receive SNRs. The 
simulation results verify our analysis and show that the rate loss per user is almost constant no matter 
where the user is located, when the channel estimators exploiting the large scale fading gains. 
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I. Introduction 

Base station (BS) cooperative transmission, which is also known as coordinated multi-point 
transmission (CoMP), is an effective way to mitigate the inter-cell interference (ICI) arisen from 
universal frequency reuse cellular systems. As a promising transmit strategy, coherent cooperative 
transmission can enhance the downlink spectrum efficiency by using multiuser (MU) multiple- 
input multiple-output (MIMO) precoding |[I1I3, when both data and channel state information 
(CSI) are gathered at a central unit (CU) via backhaul links. 

In non-cooperative systems, each BS only needs to estimate the CSI of local channels, i.e., the 
channels between the BS and the mobile stations (MSs) that are in the same cell. If the training 
sequences for the MSs in different cells are not orthogonal, the channel estimation performance 
will severely degrade due to the ICI [i3]-(3. The impact of the ICI can be mitigated by designing 
training sequences with low cross-correlation for the MSs in different cells flU, or by developing 
channel estimators exploiting the interference statistics [[51. In [[6l, the authors propose to use 
non-uniform pilot density and a DFT-based channel estimator to first separate and then subtract 
the interference signal from the estimated channel impulse response (CIR). Assuming that the 
desired channels and the interfered channels do not overlap and their interference-free initial 
estimates can be obtained through orthogonal training, the authors in [[71 propose to exploit the 
delay subspace structure to improve the estimation performance of the desired channels. 

In coherent cooperative transmission systems, the CSI of cross channels, i.e., the channels 
between the BSs and the MSs who are in different cells, needs to be estimated as well. Both the 
local and cross channel coefficients can be jointly estimated using the conventional estimators 
such as those in [8J when the training signals are orthogonal both for the MSs within a cell and 
for the MSs among the coordinated cells. However, large overhead is inevitable if orthogonal 
training signals are used for all MSs in the cooperative cell cluster. Moreover, this demands 
inter-cell signalling and protocol to coordinate the training sequences [I9l- Such a burden will 
become more noticeable when the cooperative clusters are formed in a dynamic way ifTOl . In 
[[TTl. the authors suggest to spread the orthogonal sequences from slot to slot, which may lead to 
outdated CSI at the transmitter under time- varying channels. Considering the propagation delay 
differences in multicell channels, a group of orthogonal training sequences that are robust to the 
delay are designed in [[T2l. but the number of sequences in the group is limited. 
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An inherent feature of the channels in CoMP systems is asymmetry. On one hand, the multi- 
cell downlink channels are asymmetric, which means that the average channel gains from 
different BSs to one MS are different. On the other hand, the multi-cell uplink channels are 
also asymmetric, which means that the average channel gains from MSs in different cells to 
one BS differ. Such an asymmetric channel feature is fundamental in CoMP systems, since the 
difference of the large scale fading gains cannot be compensated by an uplink or downlink 
power control mechanism. Specifically, if the MSs in different cells compensate their large scale 
fading gain differences towards one BS by power control, their receive signal energy differences 
towards other BSs will increase. This is analogous to the interference asynchrony feature, which 
cannot be dealt with by time-advanced techniques lfT3l . 

In this paper, we study the impact of the channel asymmetry on the performance of joint 
channel estimators and on the performance of downlink BS cooperative MIMO orthogonal 
frequency-division multiplexing (OFDM) systems with channel estimation errors. 

Firstly, we introduce three joint estimators requiring different channel statistics, which are the 
minimum mean square error (MMSE) estimator, a robust estimator and the least square (LS) 
estimator. We analyze the performance of these estimators when the uplink channel asymmetry 
is exploited. Our analysis shows that if the training sequences are not orthogonal among cells, 
the LS estimator will perform significantly worse than using orthogonal sequences. On the other 
hand, the MMSE and robust estimators have minor performance loss from those using orthogonal 
sequences, thanks to the large attenuation of the cross channels. 

Secondly, we analyze the impact of channel estimation errors on the performance of CoMfQ 
system using zero forcing beamforming (ZFBF) by deriving the rate loss led by the channel 
estimation errors. At the first glance, the cross channels that experience large path loss are hard 
to estimate in practice since the transmission power at the MS side is limited [|2l[T4l. which may 
degrade the downlink transmission performance. Nonetheless, our analysis shows that when the 
training sequences are orthogonal and the joint MMSE estimator is applied, the contribution of 
channel estimation errors to the rate loss is weighted by the receive SNR of the corresponding 
channel link. As a result, even though the channel estimation errors of cross channels are large, 
their impact on the rate loss is minor owing to the fact that the receive SNRs of the cross links 

'There are various transmission strategies for CoMP transmission such as coherent and non-coherent transmission. For 
simplicity, we refer the coherent BS cooperative transmission using MU MIMO precoding as CoMP transmission in the following. 
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are considerably lower than the local link. Interestingly, simulation results demonstrate that the 
rate loss of MS is nearly invariant no matter if the MS is located at the cell edge or cell center. 

The rest of the paper is organized as follows. Section II introduces the system and channel 
models. Section III and IV respectively present three joint channel estimators and analyze their 
performance. In Section V, we analyze the impact of channel estimation errors on downlink 
CoMP transmission. Simulation results are provided in Section VI to verify our analysis and to 
evaluate the system performance. The paper is concluded in Section Vn. 

Notations: Boldface upper and lower case letters X and x represent matrices and vectors, and 
standard lower case letters x denote scalars. X^, X^ and tr{X} denote the transpose, Hermitian 
conjugate transpose and the trace of X. X(i, i), X(i, :) and X(:, i) represent the element, 
the ith row and the ith column of X, respectively. ||x|| represents the two-norm of x, and diag{x} 
is a diagonal matrix with its elements. E{x} is the expectation of a random variable x. 'R{x} and 
\x\ stand for the real part and the norm of a complex scalar x. \_x\ denotes the largest integer no 
larger than a real number x and \x'\ represents the smallest integer no smaller than x. Finally, 
In denotes the identity matrix of size N, and denotes the matrix of zeros. 

II. System and Channel Models 
A. BS Cooperative Transmission System and Channel Models 

Consider a centralized CoMP system, where S BSs each equipped with Nt antennas coopera- 
tively serve M single-antenna MSs. We consider time division duplexing (TDD) systems, where 
the CSI required for MU MIMO precoding is obtained through uplink training by exploiting the 
channel reciprocity. In the uplink training phase, all MSs send training sequences and each BS 
estimates the CSI from all MSs to it. Then the BSs forward the estimated CSI to the CU via 
low latency backhaul links. The CU computes the precoding and then sends back the precoding 
vectors to each BS for downlink transmission. 

We consider frequency selective channels. The channel is assumed to be quasi-static, which 
means that the channel remains constant during the uplink training and the downlink transmission. 
The composite CIR from MS m to antenna a of BS 6 can be expressed as ^ ^ = ^ ^, 

where g^^^^ = [fi'^,b,a(0), ■ ■ ■ , gl^^b^ai^ - l)f e C^''\ am,b is the large scale fading coefficient 
including path loss and shadowing, ^ ^ = [/i^ ^ ^(0), ■ ■ ■ , f, — 1)]^ G C^^^ is the small 
scale fading channel vector, is the fading coefficient of the Ith resolvable path, which 
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is a complex Gaussian random variable with zero mean and variance a^^, and L is the number 
of resolvable paths. We assume that J2f=o = 1- 

B. Uplink Training Phase 

Except that the uplink channels are asymmetric, the signal received at one BS from MSs in 
different cells are asynchronous in CoMP systems [fT3l . Denote the propagation delay from MS 
m to BS 6 as Tm,h- We assume that the cyclic prefix in the OFDM symbol is long enough, such 
that the propagation delays turn into phase shifts in the frequency domain channels. 

Consider that all M MSs in the cooperative cluster send training sequences during the same 
uplink training duration. Denote the frequency domain training sequence of the mth MS as 
tm = [^m(O), ■ ■ ■ ,tm{K — 1)]^ G C^^^, its transmit power at each subcarrier as p^, then the 
received signal of the A;th subcarrier at antenna a of BS 6 can be expressed as 

M 
m=l 

where gl^^haW = Y^]^=q '^m,bKn,h,a{^ ~ exp(— j^//;;) denotes the composite channel fre- 
quency response (CFR) at the /cth subcarrier including both large scale fading and phase shift led 
by propagation delay, is the sampling period, K is the subcarriers number of OFDM system, 
^Lbai^) = J2f=o Knbai^) ^^Pi~j^^k) rcprescnts the small scale CFR at the A;th subcarrier, 
n{k) is the additive white Gaussian noise (AWGN) with zero mean and variance cr^. 

With the received signal, BS b estimates the composite CFR between its multiple antennas and 
all MSs, Glik) = [g(,{k),--- ,gi,,,{k'r e C^^x^*, where g;^,,(fc) = [gUM ' ' ■ ■ , oLanM^ 
is the CFR estimate between MS m and the multiple antennas of BS b. Then the BS forwards 
G((fc) to the CU, as illustrated in Fig. [Hfor a two-cell CoMP system. 

C. Downlink Transmission Phase 

The CU collects the estimated composite CFR from each BS via a low latency backhaul, and 
integrates the estimated global channel as (k) = [G{{k), ■ ■ ■ , G^^ik)] e C^'^^^^k Then, the 
CU computes the multi-cell precoding V(A;) based on &(k), and sends the downlink data and 
the corresponding precoder to each BS, as illustrated in Fig. |2l 

The downlink composite channels from multiple BSs to MS m can be expressed as g4(^) = 
[(gm,i(^))^' ■ ■ ■ ' (Sm _b(^))^]^' where gi^i,{k) G C^*^^ is the channel vector between all 



6 



antennas of BS b and MS m. Let dmik) be the data intended for MS m at the A;th subcarrier. 
For simplicity and without loss of generality, we assume that ¥.{dl^{k)dm{k)} = 1. Denote 
Vm{k) E £,BNtxi j-j^g precoding vector for MS m under CoMP transmission, and as 
the power allocated to MS m at the kth subcarrier. Then the receive signal at MS m is 

M 

y„,{k) = ^^){gUk) f^rn{k)drrXk)+ ^/^(g^ ( fc) (fc)rf, ( fc) + (fc) , (2) 

where Zm{k) is the AWGN with zero mean and variance a1 experienced at MS m. 

III. Uplink Channel Estimation for Downlink CoMP Transmission 

As shown in the downlink transmission model, the composite CFR is required for precoding 
in CoMP OFDM systems, rather than the small scale fading CFR. 

The performance of channel estimation for the composite CFR depends both on the channel 
features and on the information known a priori. In this paper, we assume that the propagation 
delays and the number of resolvable paths can be estimated perfectly. In practice, they can be 
estimated using various techniques such as those shown in [iTlfTSl . Since the number of resolvable 
paths is usually much less than the number of subcarriers in practical systems, the performance 
of channel estimation can be significantly improved by exploiting the frequency correlation of 
the channels O. When the propagation delays and the number of resolvable paths are known, 
this can simply be implemented by first estimating the composite CIR, and then obtaining the 
CFR by Fourier transformation. In the following, we only address the CIR estimation. 

To simplify our analysis, we assume that the transmit power for each MS's training sequence 
is equal. The frequency domain receive signal in ([T]) can be rewritten as a more compact form 

M M 

n,a = VP^ Y1 Tmgi,,fe,a + ^ = Vp^ T„$^,fea^,feFh^ + n 

m=l m=l 
M 

= J] X^,i,g^„,,,a + n = v^Xg*^, + n, (3) 

m=l 

where r^.^ = Ka(0), ■ ■ ■ ,n,aiK - l)f, T„ = diag{t„}, g;^;,^^ = [^^^^^(O), ■ ■ ■ , gi^i^ai^ ' 
1)]'^ denotes the composite CFR vector from MS m to the ath antenna of BS b, = 
diag{[V^„,b(0), ■ ■ ■ ,^^^b{K-l)]}, ^^,b{k) = exp(-jf ^fc), F e C^^^ is the first L columns 
of a K X K Fourier transform matrix, X^ ;, = T^^m.^F G C^^^ is the equivalent training 
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matrix of MS m by considering the known r^,?, and L, X = [Xi b, ■ ■ • , X^.b] G C^^^^^ is 
an equivalent training matrix of all MSs, g^^ = [iSi b a)^ ^ ' ' ' ^iSMba)'^V ^ C^^^^^ is the 
composite CIR vector from all MSs to antenna a of BS fe. n is the AWGN vector with zero 
mean and covariance matrix a'^Ix- 

Since each BS needs to estimate both local and cross channels for CoMP transmission, it is 
natural to estimate the CIRs from all MSs jointly, i.e., to estimate g^^. Denote the estimation 
errors of g^ as g*^„ = g^,, -g^,„ and its MSE as MSE,,„ = E{\\glJ^} = E{E11i WsUaf}- 
The MMSE estimator can be readily derived from ^ by minimizing the MSEb_a> which yields, 

gC""=(x''X+^R-)"x''r,„, (4) 

where R^ a = E{g5 ag^^} is the covariance matrix of the channels from all M MSs to the ath 
antenna of BS b. Assume that the small scale fading channels among different BS-MS links are 
uncorrected. Then Hb^a = diag{[ai^Ri,f,,a, • • • , ai/.feRA/.fe.a]}, Rm,6,a is the covariance matrix of 
the small scale fading channel vector from MS m to the ath antenna of BS b. 

Although both Rm,6,a and a^ f, vary slowly and can be estimated in practice UBl, a.^^ is 
a scalar parameter that can be estimated more accurately and requires a much lower feedback 
rate. When ^ is estimated perfectly but Rm,6,a is unknown, by assuming uniform power 
delay profile (PDP) for small scale fading channels similarly to the conventional robust channel 
estimation algorithms [fTTl . we obtain a robust channel estimator, 

g£r=(x"X + ^D-)"x%., (5) 

where D^,^ = dmg{[al^j^lL, • • • , (yli^bi^L]}- 

When we know nothing more than Tm,b and L, we can apply the LS estimator as 

gt: = (X^X)-^X%.. (6) 

IV. Performance Analysis of The Channel Estimators 

In this section, we analyze the performance of the joint channel estimators. We derive the 
MSE of the composite CIR estimates. Then we discuss the impact on the performance of the 
estimators when the training sequences are orthogonal or non-orthogonal. 
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When more than two MSs send training sequences in the uplink, it is nontrivial to obtain an 
explicit expression of the MSE of the CIR estimate. For mathematical tractability, we consider a 
simple but fundamental scenario, where B multiple-antenna BSs cooperatively serve two single- 
antenna MSs, i.e., M = 2, and the two MSs are located in two cells. 

A. MSE of Three Estimators 

We first derive the estimation error covariance matrix of the CIRs from all MSs to the ath 
antenna of the 6th BS, Rg,^ = E{g^^(g^^)^}, and then we can get the MSE of the CIR 
estimation as MSE^ b a = YlT={rn-i)L+i ^^g^ From dH), ([5]) and the covariance matrix 

for the estimators can be obtained as follows by applying the Woodbury matrix identity [ITSll . 

\ (7a) 

n / 

^j^MMSE^ (7b) 

(7c) 



where B = X^X = | _ ' ^ I ' P-,- = ^".h^m,b = F^$^_,T^T^$„,,F is the equiv- 

alent auto-correlation matrix of the training sequence for MS m, m = 1,2, Q2,i = X^^Xi ^ = 
F^$^{,T^Ti$i 5F is the equivalent cross-correlation matrix of the training sequences of MS 
1 and MS 2, and A™^-* = ^(B + ^D,-;^)-iB(Iml - RmD,-^)(B + ^D,;^)-i + R,,,B[(B + 

^Rri)-i-(B + 4D," 
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To further simplify our analysis and gain some insight into the problem, we assume uniform 
PDP of the small scale fading channels. Then the MMSE estimator degenerates to the robust 
estimator. From (TTal) . the MSE for MMSE estimator is derived as (see Appendix |A] for details) 

MSE^;^f = VmM^^ J2to ^mmse(AO, m = 1, 2, (8) 

where r/^,^ = 1/(1 + /mmse(A/) = (3 = 1/ n'=i(l + J^i) and is the Ith 

eigenvalue of < A^ < 1. 

From (TTcI) . the MSE for LS estimator can be derived as (see Appendix |B] for details) 

MSE^s,^ = ^1 ^^J^' /ls(AO, m = 1, 2, (9) 
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1-A- 



where /Ls(Ai) 

To minimize the MSE of the estimators, the matrix B should be a diagonal matrix |fT9ll . This 
requires that CI) Pmm is diagonal, which can be satisfied when T^T^ = Ik, and C2) Q2,i = 0, 
which demands the training sequences of the MSs in two cells to be orthogonal. 

The condition CI) holds when the training sequences have perfect auto-correlation. When the 
virtual carriers (VC) in practical OFDM systems are considered or the training sequences are 
not sent on all subcarriers with equi-power and equi-spaced, CI) does not hold any more. 

To ensure C2), the training sequences for the MSs in different cells should be orthogonal. This 
can be implemented in time or frequency domain, but the resources occupied by training will 
increase linearly with the number of MSs. Phase shift orthogonalization, where two sequences 
are orthogonal when their relative phase shift is larger than ^L, is known as an efficient way 
to generate training sequence with perfect cross -correlation [l8l[T9l- At most [K/L\ orthogonal 
sequences can be constructed from a sequence. 

Considering the propagation delay in multi-cell scenarios, to construct the phase shift orthogo- 
nal training sequences, the relative phase shift of the two sequences for two MSs should exceed 
^L, where L = L + Ideiay and Ideiay = [ ^'"■^^'■'' ] represents the sampled propagation delay 
difference. Then, the maximum number of orthogonal training sequences that can be constructed 
is \_K/L\. When all MSs are located in one cell, Ideiay is much smaller than L in typical outdoor 
channel^. By contrast, if the MSs are scattered in multiple cells, Ideiay will be comparable to 
the multipath delay. Consequently, few orthogonal training sequences are available for a given 
sequence length since L is large. This again leads to low spectrum efficiency. Furthermore, 
the inter-cell orthogonality demands inter-cell signalling and protocol to coordinate the training 
resourcej^, which will become a burden when the coordinated clusters are formed dynamically. 

In the following, we will analyze the performance loss led by the non-orthogonal training. To 
highlight the impact of the non-orthogonal training sequences for the MSs in different cells, we 
assume that the condition CI) holds. 

^Take the urban macro channel in systems complying Long Term Evolution (LTE) standard as an example, the multipath 
delay is usually 4 ~ 5/is. If we consider the cell radius to be 250m, then the maximum delay difference is 0.7/is, which is 
negligible compared to the multipath delay. 

^In practical cellular systems such as those complying LTE standard, partial band may be used for uplink training to increase 
the power spectrum density for cell edge MSs |20|. When there is no inter-cell coordination, the training signals may overlap 
partially in the frequency domain, which will increase the cross-correlation of the training signals 1211 . 
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B. Impact of Non-Orthogonal Training 

For comparison, we first assume that tlie training sequences of MSs in different cells are 
orthogonal, i.e., Q2,i = 0. Then the values of A;, / = 0, ■ ■ ■ , L — 1, in ([8]) and ^ are zeros, and 
we can see the MSB for estimating the local and the cross channels. In this case, both fhsi^i) 
and /mmse(Az) are equal to 1, and the MSB of MMSB and LS estimator^ are 



MSE^S" = : ^44, (10a) 



MSE^s = 1^^. (10b) 

For the LS estimator, the MSB of the composite CIR estimate, g^f,^' depends on o"^. We 
assume that the noise variance at all BSs are the same. Then the MSB of ^ ^ for 6 = 1, ■ ■ ■ ,B 
are identical, no matter they are the local or the cross channels of MS m. For the MMSB 
estimator, the MSB of g^ ^ ^ depends on ^ as well, which is the large scale fading energy of 
g^^^. If MS m is in the same cell as BS Cm, then Zm,cm,a the local composite channel for 
MS m while Zmba & 7^ are its cross composite channels. Since the local channel energy 
a^^^ is usually larger than the cross channel energy ^, h ^ Cm, we can observe from (|10al) 
that the MSB of the weak cross channels is even less than that of the strong local channels. 

At the first glance, this conclusion is inconsistent with the conventional understanding, where 
the MSB of the estimates of the cross channels should be larger than that of local channels. 
Nevertheless, this understanding is only applicable for estimating the small scale fading channels 
whose average energy is 1. To see this, we normalize the MSB of g^j„ by to obtain a 
normalized MSB (NMSB) of g^ ^ ^, which is actually the MSB for estimating the small scale 
fading channel h*„ ^ ^. The NMSB for both MMSB and LS estimators can be expressed as follows 

NMSE::^^^ = 1^^^^^, (iia) 

NMSE^s,„ = -4^-F7- (lib) 

It follows that the MSB of the estimates for the small scale fading channels with low receive 
energy is larger than that with high receive energy. 

''Again we assume unform PDF in this subsection. 
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When the training sequences of MSs in different cells are not orthogonal, then Xf ^ 0, and 
both /Ls(Ai) and /mmseIA;) exceed 1. 

From the expression of /Ls(Ai), we can see that if is close to 1 for any /, its value will be 
extremely large and the estimation performance will be severely degraded. This means that the 
LS estimator is quite sensitive to the orthogonality of the training sequences. 

In the expression of fuMSEiX), Af is weighted by /3, whose value is always less than 1. If two 
MSs are all in the same cell, is generally large, then /3 is close to 1 and /MMSE(Ai) ~ fhsi^)- 
This implies if the MSs in the same cell use non-orthogonal training sequences, the performance 
of the MMSE estimator will degrade severely. On the other hand, if two MSs are in different 
cells, since the large scale channel gains from MSs to their non-serving BSs are low that leads to 
small (3, /MMSE(Ai) will not be too large even if \f is close to 1. This indicates that the MMSE 
estimator is robust to the non-orthogonality of the training sequences in different cells, thanks 
to the severe energy attenuation of the channels from MSs to their non-serving BSs. 

Note that this conclusion holds for both MSB and NMSE since they only differ in a constant. 



C. Performance Gap between MMSE Estimator and Robust Estimator 

In wideband cellular systems, the PDP is in fact not uniform. Then the robust estimator 
will be inferior to the MMSE estimator. Nevertheless, we will show in the following analysis 
that the performance gap between the two estimators is minor when the training sequences are 
orthogonal. We will show through simulations in Section |Vl] that the same conclusion can be 
drawn when the training sequences are not orthogonal. 

When the training sequences are orthogonal, B = X^X = K\2L- Substituting B into (TTal) 
and (ITbl) . we can derive the MSE difference of the MMSE estimator and robust estimator as 

niL 

A MMSE _ 
m,b,a / J 



(robust/- ■\ "dMMSE/ 



i=(m-l)L+l 
^2 



3 



<J>;-^-^^l, (12) 



where /i = "^"u^ , crfu is the variance of the /th resolvable path of small scale fading channel. 

When ^ is large enough, /i approaches to and ^^^^^ approaches to also. On the other 
hand, when decreases, fi will increase, but Aj^^a^ will still be fairly small. That is to say, 
in the asymmetric channels, the robust estimator performs closely to the MMSE estimator. 
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V. Impact of Channel Estimation Errors on Cooperative Transmission 

In this section, we first analyze the average per MS rate loss led by the channel estimation 
errors with CoMP transmission using ZFBF. Then, we obtain a lower bound of the average 
achievable rate when the MMSE estimator is applied. 

A. CFR Estimation Errors 

Since composite CFR is required for precoding in OFDM systems, we need to transform the 
MSE of CIR estimate to the MSE of CFR estimate at each subcarrier, at (k), k = 0, - ■ ■ ,K — 
1. Denote the sum MSE of CFR at all subcarriers as E{||g^^ ;,^^||2} = E{||g^ ;,^^ - g^b.aP}- It 
can be obtained from the MSE of CIRs provided in previous sections as follows, 

nUUa\?} = E{||F„,,g^, - F^,,g^,f } = E{||F„,,g^,f } ™SE„,,,„ (13) 

where (a) comes from the fact that F^ j,Fm,f, = F^$^^$mf,F = KIl. When the training 
sequences of all MSs are orthogonal and the resolvable multipaths are uncorrelated, the MSE 
of the CFR at each subcarrier can be obtained as [fTTll 

al^^S^) = mzUaf}/K = MSE„,,,,, A; = 0, ■ ■ ■ , - 1. (14) 

From ([8]) and ^ we know that the MSE of the CIR estimates between all antennas of BS h 
and MS m are the same, which results in cr^ = crl , a = 1, ■ ■ ■ , Nt. The composite CFR 
between BS b and MS m at the A;th subcarrier can be modeled as g^^ ^{k) = g;^ ^(A;) + g;^ ^(/c), 
where gi^^ik) = [9L,b,iik), ■ ■ ■ , gi^hM^^)]" ^ sL,bik) is the estimation of gi^^ik) and g^^^ik) is 
the estimation error vector whose covariance is al^^{k)Ij^^. Since the transmission procedures 
of all subcarriers are same, the index of subcarrier is omitted in the following for brevity. 

B. Impact of Channel Estimation Errors on CoMP Transmission 

When the global channel vectors are reconstructed at the CU from the estimates provided by 
all coordinated BSs, a multicell ZFBF is computed as follows 

(15) 

Then the beamforming vector of all cooperative BSs for MS m is obtained by normalizing the 
mth column of V as v„ = V(:, m)/|| V(:, m) || and v„ = [(vi,„)^, ■ ■ ■ , (v^^^)^]^, where 
\'b,m e C^*^^ is the precoder vector of BS b for MS m. 
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In order to derive a closed-form expression of the per MS rate loss led by the channel 
estimation errors, we assume that the number of MSs cooperatively served by B BSs is BNt, 
which indicates full multiplexing CoMP-MU transmission as in [22]. We further assume that the 
power allocated to all MSs are identical, which is denoted as p^. 

The average rate of MS m achieved by CSI estimate-based ZFBF is obtained from ^ as 



i?„ = E {log2(l + SINR^)} = E <; log2 I 1 + /"Jg^^-I' ) K . (16) 




The average rate of MS m achieved by perfect CSI-based ZFBF is given by 

oldeal Ti? fi , P \„H Ideal 12^1 

Rm = E |l0g2 yl + -^\Zm^m I jji (17) 

where vj^^^' = [(v^';'^^')^, ■ ■ ■ , (v^_^^')^]^ is the perfect CSI-based ZFBF vector of all BSs for 
MS m, which is chosen to be orthogonal to for j = 1, ■ ■ ■ , M and j ^ m. 

Theorem 1 The rate loss of MS m of the CoMP transmission using ZFBF led by the channel 
estimation errors can be upper bounded by 

Ai?™ = i?^^^' -Rn.< log2 (^1 + = Ai?UB, (18) 

where = Z]j=i is th^ average interference power experienced by MS m, 

§m = [(gm.i)^! ■ ■ ■ ! (Smfi)^]^ estimation error vector of the global channel vector g^- 

The derivation is similar to that in [|22ll . Due to the lack of space, we omit the proof of the 
Theorem. From (fTSi) . the achievable rate of MS m when estimated CSI are used for CoMP 
transmission can be lower bounded by 

> R':^'^' - Ai?UB. (19) 

To gain further insight into the rate loss, we assume that the channel estimation errors g^ 
are independent of the precoder vectors for j = 1, ■ ■ ■ , M and j ^ m. This assumption is 
satisfied when the MMSE estimator is applied. When MMSE estimator is used, since the channel 
estimation errors are independent of the channel estimates, and the precoders are functions of the 
channel estimates, the channel estimation errors and the precoders are mutually independent. In 
Section |Vll we will examine the impact of estimation errors led by the LS and robust estimators 
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on CoMP transmission through simulations. 

As stated in Section IV-Al the MSEs of the CFRs between MS m and all antennas of BS b are 
identical, then the rate loss upper bound of MS m can be further derived as follows by taking 
expectation over the channel estimation errors (see Appendix O for details) 



AR 



UB 



M -B 2 
j=ljjtm 6=1 




V6,' 



(20) 



Remark 1 To connect with the conventional understanding of the impact of the channel esti- 
mation errors, we show the contribution of channel estimation errors of the small scale fading 
channels, NMSE^ to the rate loss in (|20l ), which is weighted by the downlink receive SNR 
of the link SNR^ ^. For the local channels of MS m, i.e., b = Cm, NMSE^ small and 
its contribution to the rate loss will be minor. For the cross channels of MS m, the estimation 
errors NMSEm,6 for b ^ Cm will be large. However, because the receive SNRs of the cross 
links are considerably low, the impact of the channel estimation errors of cross channels will be 
significantly alleviated. This is true especially for cell center MSs. 



Substituting the NMSE^.t of the MMSE estimator under orthogonal training shown in (lllal) 
into (|20|) . we obtain the rate loss upper bound as follows 
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_i .-i™ J, 1 -L 



V6,1 



< log2 



where (a) is obtained because 



< 1, and YX=i I|v6,ii 



(21) 



^ ' alp"" K 



1 as described in (1151) . p 



and are respectively the power transmitted to each MS and that transmitted by each MS, a1 
and cr^ are the noise variances at the BS and MS. When p'^ and are fixed, the upper bound 
of the rate loss will not depend on the large scale fading gains of both local and cross channels. 

According to (fT9l ) and (|2TI) . we can obtain the lower bound of the average rate achieved by 
MS m under MMSE estimator and orthogonal training as 



Rm > R: 



Ideal 
m 



logs 



1 + (M 



cr2 p« K 



R, 



LB 



(22) 
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VI. Simulation Results 

In this section, we compare the performance of different channel estimators and evaluate their 
impact on the performance of downlink CoMP system. 

A cooperative cluster of two cells is considered (B = 2). The cell radius r = 250 m. Each 
BS has four omnidirectional antennas (Nt = 4), serving two single- antenna MSs. The downlink 
transmit power for each MS p'^ is 5 dB larger than the uplink transmit power of each MS p". The 
maximum delay spread of the channel r = Afis, and the channel is implemented as a tapped-delay 
line with Rayleigh fading coefficients and an exponential DPD with the attenuation factor being 
1.4. K = 128. The system bandwidth is 5 = 5 MHz, the sampling period is = 1/5 = 0.2/iS, 
thus L = t/Ts = 20. Due to the lack of the space, we only present the performance of an 
OFDM system without VC and the training sequences are sent on full band. Extensive results 
show that the same conclusions can be drawn to the OFDM system with VC and the training 
sequences of MSs are sent on partial band. 

The training sequences are constructed from Constant Amplitude Zero Autocorrelation Code 
(CAZAC) (231 as t{k) = e'^ ^zc , k = O,-- - ,K - 1, [20], where A^zc = 127, Uk = 
mod(/c, Nzc)- The training sequences for MSs in the same cell are orthogonal by cyclic shifting. 
The training sequences of MSs in different cells can be orthogonal or non-orthogonal. For 
orthogonal training, the training sequences for the four MSs in the two cells are constructed 
from the cyclic shift of the CAZAC with the same value of c. For non-orthogonal training, the 
training sequences of MSs in two cells use different values of c. The values of c for two cells 
are set to ci = 1 and C2 = 7, considering that their cross correlation is moderate. 

A. NMSEs of Dijf event Estimators 

To show the impact of MSs' positions on the estimation errors for small scale fading channels 
under orthogonal and non-orthogonal training for multiple MSs, we let the four MSs in the two 
cells be symmetrically located, as shown in Fig. [3l Then the channel estimation performance of 
all MSs are the same. We take the performance of one MS as an example to analyze. 

In Fig. m the NMSEs versus local uplink receive SNR of three estimators for both local and 
cross channels are shown. When the training sequences are not orthogonal, the performance of 
the LS estimator degrades severely. By contrast, the performance gap of the MMSE estimator 
under orthogonal and non- orthogonal training is minor. Comparing the NMSE of the robust and 
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MMSE estimators, we observe that the performance loss of the robust estimator from the MMSE 
estimator is small. These results agree well with our previous analysis. Again, we should note 
that the impact of the non-orthogonal training sequences on the performance of estimating the 
small scale fading channels is the same as that of estimating the composite channels, since the 
MSE and NMSE only differ in a constant. 

B. Downlink Average Rate with Different Channel Estimators 

1) Positions ofMSs are Fixed: Now we verify the analysis in Section V, where the positions 
of the MSs in two cells are the same with before. We simulate the downlink average throughput 
of each MS, which is averaged over 1000 realizations of small scale fading channels. 

We first evaluate the tightness of the rate lower bound derived in Section |Vl The performance 
with MMSE estimator under orthogonal training is taken as an example. The per MS rate under 
the pefect ZFBF and the CSI estimate-based ZFBF are shown in Fig. [51 together with the lower 
bound of the achievable rate derived in (|22l). It shows that the derived lower bound is close to 
the rate obtained by simulation. 

We then compare the impact of different estimators on the performance of CoMP transmis- 
sion under both orthogonal and non-orthogonal training, which is shown in Fig. [6l where the 
throughputs under Non-CoMP transmission are also present as a reference. When the training 
sequences are orthogonal, the per MS throughputs under different estimators are almost the same. 
The performance gap from the perfect CSI-based ZFBF does not change no matter the MS is 
located at the cell edge or the cell center. When the training sequences of MSs in different 
cells are not orthogonal, the performance degradation when using both the MMSE and the 
robust estimators is minor. On contrary, the performance when using the LS estimator severely 
degrades, which is even worse than the Non-CoMP transmission. It indicates that the estimator 
can perform fairly well only when the large scale fading information is employed. 

2) Positions of MSs are Randomized: Finally we simulate the case in which the locations of 
two MSs are randomly distributed in each cell. We use the average throughput per MS and the 
cell edge MS throughput as the performance metric. The results are shown in Fig. |71 which is 
obtained from 1000 random drops. We can see that the same conclusion can be drawn as the 
case where the positions of MSs are fixed. 
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VII. Conclusion 

In this paper, we have studied the impact of uplink channel asymmetry on the performance of 
channel estimation and the impact of downlink channel asymmetry on the performance of BS 
cooperative transmission with channel estimation errors. We have analyzed three joint channel 
estimators, the MMSE estimator, a robust estimator and the LS estimator. Our analysis showed 
that if the training sequences of MSs in different cells are not orthogonal, the performance 
of estimating cross channels with the LS estimator severely degrades. By contrast, the MMSE 
estimator is robust to the non-orthogonal training, and due to the uplink channel asymmetry the 
robust estimator has minor performance loss. By analyzing the impact of channel estimation 
errors on the cooperative ZFBF transmission, we showed that due to the downlink channel 
asymmetry, the contribution of the channel estimation errors to the rate loss is weighted by the 
receive SNRs of the corresponding links. As a result, despite that the estimation errors of the small 
scale fading channels of the cross links are large, their impact on the rate loss is minor. When 
the joint channel estimators exploit the large scale fading gains, CoMP transmission performs 
fairly well, even without inter-cell orthogonal training. This improves the spectral efficiency and 
simplify the inter-cell signalling required to coordinate the training resources. 

Appendix A 
Derivation of the MSE of MMSE estimator 

Assume uniform PDP for small scale fading channels, i.e., R^^a = diag{[aj f^jlL,(y2^f^jiL]}- 
Substituting the expression of B to (TTal) and applying the formula of block matrix inversion, the 
covariance matrix of the estimation errors becomes 



a' 



Rrr' = :^^l ^ . ..^^.^ I' (23) 



where N = (^I, - v,,2^) \ V,,m = 



2 

b.m^ K 

The MSE of the MMSE estimator for CIR from the MS 1 to antenna a of BS 6 is 



MSE--- = ^tr{N} = ,,,^tr | (l, - ,,,,,, 



Q21Q 



21 



L-1 

(a) 
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where (a) follows from applying the eigenvalue decomposition as ^'^^^^ = UAU^, U is an 
unitary matrix and A = diag{[AQ, ■ ■ ■ , A|_J}. 

The MSE of MMSE estimator for CIR from the MS 2 can be derived as 



^MMSE "^n +^ f_2 Q21-|vtQ21 



2 J. I f ^ T Q2lQ2l\ Q2I Q '"^ 



n 2 + J / T ^21^21 \ ^21 ^21 ( , "n + fT 1 



0"„ 



2 L — 1 ^9 9 9 L — 1 

„2 



T.-r^ + ^^M^ = ^^2^E^^Tl^f^ — • (25) 



Comparing (|24|) and ([25]) . we can see that the only difference is the factor r/bm- The channel 
that experiences large attenuation exhibits small estimate errors. 

Appendix B 
Derivation of the MSE of LS estimator 

Substituting the expression of B into (|7c1) and using the formula of 2 x 2 block matrix inversion, 
we can derive the covariance matrix of the estimator errors as 



where M= (^I^- Q^)"'. 

The MSE of the LS estimator for the CIR from MS 1 to antenna a of BS 6 can be derived as 



L-l 



MSEg, = ^^-{M} = ^tr { (I, - UAU-)-^} = (^7) 
The MSE of LS estimator for CIR from the MS 2 can be derived in the same way as 



pu^ [ K p'^K 

. \ fr Q2lQ2l\ Q2I Q 



L-l 

(28) 

I — AT 

1=0 



Comparing dJV]) and dlH]), we find that the MSEs of the LS estimator for the CIRs from 
different MSs are identical. 
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Appendix C 

Rate loss after taken expectation over channel estimation errors 

The average interference power experienced by MS m in (fTSl ) can be derived by taking 
expectation with respect to the estimation errors, i.e., 

2" 



E{Im} = E.f 




M 



&m 



m.,b) ^b,j 



E E 



(29) 



j=ljjtm b=l 

where (a) follows because the channel estimation errors of the channels from multiple BSs 
to one MS are assumed to be uncorrelated and their expectations are zero, (b) is derived by 
considering that the covariance of estimation errors for the channels from multiple antennas of 
one BS are the same, i.e., E{gl0l ,,)"} = <^ ,IjVf 

Substituting (|29l) into (fTSi) . we can get the upper bound of the rate loss under CoMP trans- 
mission as 

rl M B 



AR. 



UB 



log2 



1 + 



p 



■ E E 

M 



1 1 \ ^ '^m,bP ^em,b II ||2 

1+ ^2 .,2 11^^.-" 



0"f 



a: 



m,b 



(30) 



j=lj=^m 6=1 

References 

[1] M. K. Karakayali, G. J. Foschini, and R. A. Valenzuela, "Network coordination for spectrally efficient communications in 

cellular systems," IEEE Wireless Commun. Mag., vol. 13, no. 4, pp. 56-61, Aug. 2006. 
[2] A. ToUi, M. Codreanu, and M. Juntti, "Cooperative MIMO-OFDM cellular system with soft handover between distributed 

base station antennas," IEEE Trans. Wireless Commun., vol. 7, no. 4, pp. 1428-1440, Apr. 2008. 
[3] J. Jose, A. Ashikhmin, T. Marzetta, and S. Vishwanath, "Pilot contamination problem in multi-cell TDD systems," in Proc. 

IEEE Int. Symp. Information Theory (ISIT), Seoul, Korea, June 2009, pp. 2184-2188. 
[4] G. Kang, P. Hasselbach, Y. Yang, P. Zhang, and A. Klein, "Pilot design for inter-cell interference mitigation in MIMO 

OFDM system," IEEE Commun. Lett., vol. 11, no. 3, pp. 237-239, Mar. 2007. 



20 



[5] D. Katselis, E. Kofidis, and S. Theodoridis, "On training optimization for estimation of correlated MIMO channels in the 

presence of multiuser interference," IEEE Trans. Signal Processing, vol. 56, no. 10, pp. 4892^904, Oct. 2008. 
[6] S. Lee, K. Kwak, J. Kim, and D. Hong, "Channel estimation approach with variable pilot density to mitigate interference 

over time-selective cellular OFDM systems," IEEE Trans. Wireless Commun., vol. 7, no. 7, pp. 2694—2704, July 2008. 
[7] M. R. Raghavendra, S. Bhashyam, and K. Giridhar, "Interference rejection for parametric channel estimation in reuse-1 

cellular OFDM systems," IEEE Trans. Veh. TechnoL, vol. 58, no. 8, pp. 4342-4352, Oct. 2009. 
[8] Y. Li, "Simplified channel estimation for OFDM systems with multiple trasnmit antennas," IEEE Trans. Commun., vol. 1, 

no. 1, pp. 67-75, Jan. 2002. 

[9] H. S. Kim, S. H. Lee, and Y. H. Lee, "Channel sounding for multi-sector cooperative beamforming in TDD-OFDM wireless 

systems," in Proc. IEEE Int. Conf. Commun. (ICC), Cape Town, South Africa, May 2010. 
[10] A. Papadogiannis, D. Gesbert, and E. Hardouin, "A dynamic clustering approach in wireless networks with multi-cell 

cooperative processing," in Proc. IEEE Int. Conf. Commun. (ICC), Beijing, China, May 2008, pp. 4033-4037. 
[11] M. S. L. Thiele, S. SchiffermuUer, V. Jungnickel, and W. Zirwas, "Multi-cell channel estimation using virtual pilots," in 

Proc. IEEE Veh. Technol. Conf. Spring (VTC-S), Singapore, May 2008, pp. 1211-1215. 
[12] T. Kwon, H. Song, and D. Hong, "Robust channel estimation in multicell OFDM(A) downlink systems with propagation 

delay," in Proc. IEEE Veh. Technol. Conf. Spring (VTC-S), Dublin, Ireland, May 2007, pp. 1450-1454. 
[13] H. Zhang, N. B. Mehta, A. F. Molisch, J. Zhang, and H. Dai, "Asynchronous interference mitigation in cooperative base 

station systems," IEEE Trans. Wireless Commun., vol. 7, no. 1, pp. 155-165, Jan. 2008. 
[14] S. Han, C. Yang, M. Bengtsson, and A. I. Perez-Neira, "Channel norm based user scheduling in coordinated multi-point 

systems," in Proc. IEEE Glob. Telecom. Conf. (GloheCom), Nov. 2009. 
[15] M. R. Raghavendra and K. Giridhar, "Improving channel estimation in OFDM systems for sparse multipath channels," 

IEEE Signal Processing Lett., vol. 12, no. 1, pp. 52-55, Jan. 2005. 
[16] K. C. Hung and D. W. Lin, "Pilot-based LMMSE channel estimation for OFDM systems with power-delay profile 

approximation," IEEE Trans. Veh. Technol., vol. 59, no. 1, pp. 150-159, Jan. 2009. 
[17] Y. Li, L. J. Cimini, and N. R. SoUenberger, "Robust channel estimation for OFDM systems with rapid dispersive fading 

channels," IEEE Trans. Commun., vol. 46, no. 7, pp. 902-915, July 1998. 
[18] G. H. Golub and C. F. V. Loan, Matrix Computations. Baltimore MD: Johns Hopkins University Press, 1996. 
[19] 1. Barhumi, G. Leus, and M. Moonen, "Optimal training design for MIMO OFDM systems in mobile wireless channels," 

IEEE Trans. Signal Processing, vol. 51, no. 6, pp. 1615-1624, June 2003. 
[20] 3GPP Long Term Evolution (LTE), "Physical channels and modulation," TSG RAN TR 36.211 v8.4.0, Sept. 2008. 
[21] Y. Ogawa, T. Takata, T. Iwai, D. Imamura, K. Hiramatsu, and K. Miya, "Pilot signal generation scheme using frequency 

dependent cyclic shift sequence for inter-cell interference mitigation," in Proc. IEEE Radio and Wireless Symposium (RWS), 

San Diego, CA, Jan. 2009, pp. 421^24. 
[22] G. Caire, N. Jindal, M. Kobayashi, and N. Ravindran, "Multiuser MIMO achievable rates with downhnk training and 

channel state feedback," IEEE Trans. Inform. Theory, vol. 56, no. 6, pp. 2845-2866, June 2010. 
[23] D. C. Chu, "Polyphase codes with good periodic correlation properties," IEEE Trans. Inform. Theory, vol. 18, no. 4, pp. 

531-532, July 1972. 



21 




Fig. 1. Uplink training procedure of a two-cell CoMP system, where the solid lines denote the local channels and the dash 
lines represent the cross channels. The frequency domain channels are shown and the index of subcarrier is omitted for brevity. 




Fig. 2. Downlink transmission procedure of a two-cell CoMP system. 




Fig. 3. MSs' positions in the simulated CoMP system. The locations of two MSs in the same cell are symmetrical to the 
line connecting the two BSs, and the MSs in different cells are symmetric to the cell edge, ds is fixed to be ^. All MSs 
move from the cell edge to the cell center simultaneously, then their local receive SNRs all increase. Given the value of di, 
we can get the value of d2 and vise versa. Assume that the downlink receive SNR of the cell edge MS, SNRedge, is lOdB. 
Consider the path loss factor e as 3.76, then the receive SNR of a MS from a BS with a distance d can be computed as 
SNR(d) = SNRedge + £l01ogio(5). Similaiiy, when SNR(d) is given, we can get d. 
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Local SNR^'-(dB) 



Local SNR^'-(dB) 



Fig. 4. NMSEs of different estimators for both local and cross channels versus the receive SNR of the local channel. The NMSEs 
are obtained by averaging over 1000 realizations of small scale fading channels. The X-axis is defined as Local SNR^"^ = 
, where afocai is the large scale fading energy of the local channel for the MS. When the MSs are at the cell edge, 
Local SNR"^ = 5 dB. The Y-axis is the NMSE of channel estimators, which reflects the estimation performance of the small 
scale fading channels. When Local SNR'^'" increases, the large scale fading gains of the cross channels decrease, which leads 
to large NMSE of the cross channels. For the local channels, the NMSE of the three estimators are overlapped under orthogonal 
training (shown as "-0" in the legend). For the cross channels, the NMSE of the MMSE estimator under non-orthogonal training 
(shown as "-NO" in the legend) is overlapped with that of the robust estimator under orthogonal training. 




Fig. 5. Achievable rate and its lower bound of a MS when the MMSE estimator and orthogonal training are considered. The 
per MS rate of CoMP transmission with perfect CSI is also provided for reference, which is shown as "ideal CSF' in the legend. 
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Fig. 6. Achievable rate of a MS when different estimators are used with both orthogonal and non-orthogonal training. For 
the Non-CoMP transmission with estimated CSI, the CSI for downlink preceding is estimated by the conventional single user 
MMSE estimator |8|. The performance when the robust estimator and the LS estimator under orthogonal training are applied 
overlap with that when MMSE estimator under non-orthogonal training are used. The meaning of the legends is the same as 
previous figures. 




Fig. 7. Cell average and cell edge MS throughput when different channel estimators are used. The cell edge MS throughput is 
defined as the 5% point of the cumulative distribution function of the MS throughput. As a baseline, the results of Non-CoMP 
transmission are also provided. The meaning of the legends is the same as previous figures. 



